Dynamics and stress in gravity driven granular flow 
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We study, using simulations, the steady-state flow of dry sand driven by gravity in two-dimensions. 
An investigation of the microscopic grain dynamics reveals that grains remain separated but with a 
power-law distribution of distances and times between collisions. While there are large random grain 
velocities, many of these fluctuations are correlated across the system and local rearrangements are 
very slow. Stresses in the system are almost entirely transfered by collisions and the structure of 
the stress tensor comes almost entirely from a bias in the directions in which collisions occur. 
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Despite being the raison d'etre for one of the oldest 
scientific instruments, an understanding of the dynam- 
ics of sand grains in the hourglass remains a mystery 
which still leaves much to be explored. Recent research 
into granular material has followed a number of differ- 
ent tracks [0. The initial state of a diffuse granular gas 
with velocity fluctuations (about some mean flow veloc- 
ity (v)) drawn from a Maxwcll-Boltzmann distribution 
with characteristic "granular" temperature T, can be de- 
scribed with fluid-mechanical equations using the ideas 
and machinery of kinetic theory Yet, left to them- 

selves, such systems are unstable to the formation of high 
density clusters [|]||]. As the system cools, due to the ki- 
netic energy of the grains being dissipated in the inelastic 
collisions, correlations in the density and spatial velocity 
build up and hydrodynamic descriptions based on molec- 
ular chaos (mean field assumption) break down If 
inelastic collapse to a static state is to be avoided, the 
system must be constantly supplied with energy. 

Here, we consider inelastic hard spheres under the in- 
fluence of gravity in a cylindrical container (chute flow) 
with a finite probability of reflection p, at the bottom 
(equivalent to a sieve in an experiment). The dynamics 
conserve momentum and the energy loss in a collision is 
proportional to 1 — /x 2 where \i is the coefficient of normal 
restitution. [The grain dynamics are the same as those 
used in 10 and others]. Due to the flow restriction 
imposed by the sieve at the bottom, the density is very 
high (around 90% of that of a closed packed configura- 
tion) . The mean velocity profiles are very similar to those 
found in quasi-static flows 1 1 12] , a fairly flat region in 



the middle with exponential boundary layers, or shear 
zones at the edges. However, as there is no friction in 
the system and the grains remain separated, a different 
explanation is required here. We believe the ability of 
the system to withstand shear in this regime is due to 
the very long time scales related to the slow glass-like 
rearrangement of grains. Diffusion is extremely limited, 
something that is reflected in the finding that both the 
distribution of distances and times between collisions is 
a power law. In addition, we find that the stress tensor is 
not simply related to gradients in the density or velocity 



fields, as suggested by hydrodynamic-type descriptions of 
more diffuse granular flows |^,|l3), but obtains its struc- 
ture primarily from an anisotropic distribution of colli- 
sion impact directions. The structure of the stress tensor 
is thus found to be similar to recently proposed models 
of stress propagation in static sand in which the direc- 
tion of the characteristics for stress propagation is fixed 
|l4| . This gives insight into how such models of the static 
case may be obtained from the dynamics which were in- 
volved in the creation of static configurations of grains. 
Finally, we find large velocity fluctuations in the bulk re- 
gions which are free from shear, in agreement with recent 
experiments Jig ]. This apparent contradiction with the 
more accepted view that velocity fluctuations should be 
zero in the absence of shear [ 2|[[3[ can be resolved when 
one takes into account heat flow. It is found that veloc- 
ity fluctuations are generated in the shear zones and are 
then conducted to the interior region. 

Figure |l|(a) shows a snapshot of one of our 2-d simu- 
lations. When a grain hits the bottom of the chute, it is 
reflected with probability p and escapes with probability 
1 — p. Whenever a grain escapes from the bottom a new 
grain is placed at the top. The grains are polydisperse 
with the diameter of the ith grain being di = 2 + e f z<5j 
where the 5i are independently distributed Gaussian ran- 
dom variables with unit standard deviation, and td typi- 
cally around 0.1 — 0.2. This polydispersity breaks up the 
long range crystalline order which would otherwise be 
present (crystalline order still persists over a distance of 
3-5 grains). When two grains collide, the velocities after 
collision r[ and r' 2 , expressed in terms of the velocities 
before collision, ri and r 2 , are 
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where q = (r 2 — ri)/|r2 — i"i|, and /i is the coefficient of 
restitution (typically /i = 0.95 or fi — 0.9 in our simula- 
tions). While energy is lost to the internal energy of indi- 
vidual grains, momentum is conserved. For < p < 1, a 
steady state situation is achieved where a dense column 
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of granular material flows down the chute with a constant 
(in time) average velocity. Within this column, grains 




FIG. 1. (a) Sand-flow in a 2- dimensional cylindrical geom- 
etry. Gravity is in y-direction. The color codes the i/-velocity 
varying from yellow for fast to blue for slow. Grains are circu- 
lar with 10% polydispersity, coefficient of restitution /i = 0.95, 
and the probability of reflection at the bottom is p = 0.5. 
Profiles across the chute at y = 100 of: (b) Average vertical 
velocity (v y ); (c) Root mean square velocity fluctuations for 
the x (solid) and y-components (dashed); (d) Average density 
p as a function of x. 

remain separated, flying from collision to collision with 
fairly large velocities. The grains advance between col- 
lisions following Newton's equations of motion. Around 
10 9 collisions were tracked (~ 10 6 per grain) with an 
additional 10 s collisions discarded at the beginning to 
ensure we had reached steady state conditions. During 
this time the grains travel through the system about 20 
times. The width of the chute (50 or 100 radii) is compa- 
rable or larger than the smallest dimension of situations 
studied experimentally in 3-dimensions in |fl5f . 

If, as in the experiments, the pressure is to be inde- 
pendent of the height of the column the side walls must 
in some way be able to support the weight or absorb mo- 
mentum in the vertical direction. In order to do this we 
make collisions with the side walls inelastic in the tan- 
gential direction so that after colliding with a side wall, 
the velocity r' in terms of the velocity before the colli- 
sion is r' = {—v x ,pi y ) (There is nothing special about 
the choice of \x for the factor by which r y is decreased). 
One can think of this as a very simple model of "soft" 
walls (softer than the grains, e.g. Plexiglas walls and 
brass balls) or an even simpler model of a rough surface. 
As we shall see, the boundary layers are exponentially 
small so that the precise boundary conditions should not 



affect grain motion in the interior. 

Figure 0(b) plots the average vertical velocity in a 
cross-section of the chute. The velocity is flat in the 
central region, implying an ability to withstand finite 
shear stresses, with exponential boundary layers, or shear 
zones. These profiles are very similar to those seen in 
slow and quasi-static flows [|l3|,[ll],[L2| . However, the ex- 
planation for quasi-static flows involves overcoming the 
yield stress which arises from static friction ]l2| ). Here 
there is no friction between the grains and the grains re- 
main separated (and undergo fairly large velocity fluctua- 
tions, Figure ||(c)). The ability to withstand shear is less 
surprising when you consider correlations in the veloc- 
ity fluctuations. Figure [|(a) shows the velocity-velocity 
correlation function for the x component of velocity. As 
can be seen, the velocity correlation function scales with 
system size. As we shall see below, this correlated mo- 
tion allows for the coherent transfer of momentum across 
the system from collision to collision, thereby resulting 
in some of the solid-like properties. 

Before examining the process of stress transfer, it is 
worthwhile to look at the underlying microscopic dynam- 
ics, and how these differ from those in a normal fluid or 
solid material. In diffuse regions, the mean free path can 
be a very useful quantity to measure JlOp . Figure ||(b) 
shows the (unnormalized) probability N(t) that a grain 
survives a time t without suffering a collision, measured 
near the center of the dense column. As can be readily 
seen, N(t) has a power law dependence on r: 

N(r) ~ r- a (2) 

where a — 2.75 ± 0.05. There is a short-time cut-off 
where N(t) — > constant as r — > 0. In addition, the 
distance traveled between collisions I is also a power law 
with the same exponent. This is not that surprising since 
the moments of the velocity distributions appear to be 
convergent and thus r and I are related by some typi- 
cal velocity. The power law leads to a typical mean free 
path comparable to the short distance cutoff given by 
the inter-particle spacing. This suggests that particles 
are tightly confined by their neighbors. 

To examine this, it is instructive to look at the ob- 
served diffusion of a grain as it travels through the sys- 
tem. Figure ||(c) shows the average distance traveled in 
the ir-direction by a single particle (solid line) in a time 
t and the distance between two particles (dashed line) 
which started out as nearest neighbors near the top of 
the column. One sees that diffusion is very limited and 
that the two particles typically remain nearest neighbors 
throughout their lifetime in the column. Even the single 
particle motion is fairly limited moving only of order 1 
radius during it's lifetime in the column. The single par- 
ticle motion is not diffusive (diffusive motion would have 
(\x(t) — £c(0) |) ~ t 1 / 2 ) but appears to be approaching this 
in the long time limit. 

Momentum conservation requires that 

dt(pvi) + dk(pviV k ) = dk(Ji k + pgi (3) 
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where an, is the stress tensor, and all Vk refer to first 
moments of the velocity distribution. In our steady state 
situation the left hand side is zero. There are two types 
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FIG. 2. (a) Correlation function for the x component of ve- 
locity shown as a function of the scaled variable (x — xo)/W 
where W is the width of the system (W = 50 for the 
triangles and 100 for the squares) and measured from 
(xo,yo) = (0, 133), a point at the middle of the column, (b) 
Unnormalized probability N(t) that a ball survives a time r 
without suffering a collision, as measured near the center of 
the column. The straight line indicates the power law r -2 ' 75 . 
(c) For a grain starting at r(0) = (0, 180), the average dis- 
tance (solid line) a grain moves in the x direction as a func- 
tion of time; distance between two particles which started out 
as nearest neighbors (dashed) and; power law i 1 ' 2 expected 
for normal diffusion (dot-dashed). In (b) and (c) data is from 
a system of width 50, p = 0.5, and [i — 0.9. (d) When two 
grains collide momentum is transferred in the x direction 
across the intervening plane marked by the dashed line. 

of contributions to the stress tensor. The first contri- 
bution comes from the diffusion of particles (By this we 
do not mean macroscopic convection, the dk(pViVk) term, 
but the microscopic particle diffusion, present even in the 
case of no macroscopic flow). This contribution turns out 
to be negligible as particles do not diffuse around very 
much (Fig. |](c)). Unlike the often studied granular gas 
case, in the situation we study here the stress is domi- 
nated by momentum transfer during collisions. 

Figure ^(d) shows two grains at the moment of col- 
lision. The collision results in momentum being trans- 
ferred in the x' direction across the intervening plane 
marked by the dashed line in the Figure. This gives a 
contribution to the stress tensor, measured in the (x', y') 
coordinate system of 
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jrain 2. Converting to xy-coordinates this is 
_ / (p c • x)(x • x') (p c -y)(x-x') 
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The net contribution to the stress tensor is found by av- 
eraging over time and dividing by the area over which 
the average is performed. Figure ||(a) and (b) shows a yy 
and a xy for one particular case. Figure ||(c) shows how, 
within statistical errors, dkaik + pg% = thus satisfy- 
ing Eq. (||) (this relationship should be interpreted in the 
probabilistic sense: It is true "on average" ) . Interest- 
ingly, most of the weight is supported by a shear stress, 
via d x a xy rather than the pressure gradient d y a yy w 
0.014. Effectively, the column behaves like a solid block 
sliding down a tube with friction at the walls supporting 
the force due to gravity. However, the system has solid 
properties only in the sense of a glass, as rearrangements 
do occur, albeit on longer time scales. 

Given that things are not isotropic, it is important to 
see where the structure of the stress tensor comes from. 
It can arise from both gradients in velocities or from an 
anisotropic distribution of collision directions. If at every 
collision one draws a vector from the center of one grain 
to the center of the other grain and then averages, one 
finds that there is a definite bias in the collision directions 
(Figure [|(b), right axes). This is consistent with the colli- 
sion "chains" observed in the simulations and previously 
by other workers (|||. This bias in the collision direc- 
tions implies a preferred direction to the collision chains. 
To clarify, we put the momentum transfer from Eq.(Q) 
into the stress tensor from Eq.(||) to get 

a v=t XI -^(l+/i)(fi-f 2 )'qx 

collisions 

(q-x) 2 (q-x)(q-y) 
(q-y)(q-x) (q • y) 2 



? (l + /i)/ c ((f 1 -r 2 ).q) x 

((q-x) 2 ) <(q.x)(q.y)) 
((q.y)(q.x)) ((q-y) 2 ) 



(6) 



where p c is the momentum transfered from grain 1 to 



where the sum is over collisions in a (long) time interval t 
and f c is the collision frequency. In writing this we have 
made an assumption, supported by the data, that the fac- 
tor (fi — 1*2) • q can be separated from the factors in the 
matrix when computing averages. In fact, the factors in 
front of the matrix appear to be nearly a constant, hence 
the structure of the stress tensor comes almost entirely 
from the bias in the directions in which the collisions oc- 
cur (cf. Fig. ||(b)). Note that there is a striking similarity 
between the structure of the stress tensor in our model 
and that of a static model recently proposed in Ref. Q . 
In our case, the dynamics of flow automatically drive the 
system to a "jamming" configuration where force chains 
(as given by the preferred directions of collision) form and 
the system acquires an ability to withstand shear stress. 
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Finally, we examine the energy balance. As we saw 
in Figure 0(c), there are fairly large velocity fluctuations 
throughout the flow. This is in agreement with the ex- 
periments of Menon and Durian jig] , but in apparent 
contrast to hydrodynamic descriptions of granular flows 
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FIG. 3. (a) Diagonal components of the stress tensor a yy 
computed from momentum transfer during collisions, (b) 
Off-Diagonal component of the stress tensor, a yx (solid line 
and left axes) and correlation between x and y components 
of a vector q drawn at every collision from the center of one 
grain to the center of the other grain (dashed line and right 
axes), (c) Balancing of body force pg (solid line) and deriva- 
tives of stress tensor (d x a xy +d y a yy ) (dotted line). The larger 
fluctuations in pg near the edge are due to columnar ordering 
induced by the straight side walls. Statistical error is compa- 
rable to the noise level in (d x a xy + d y a yy ). Data is from the 
system of width 50, p = 0.5, and p = 0.9. (d) Heat flow in 
^-direction Q x (solid) and —pfc\d x (v 2 ) (dotted). 



13 which suggest that velocity fluctuations about the 
mean flow velocity should be zero in the absence of shear. 
A resolution to this apparent conflict may be found by 
examining the energy flow. For our steady state situa- 
tions, energy conservation requires 



= -V • Q + pg 
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From Eq. (jl]), the (kinetic) energy lost in a single col- 
lision is 8E = — j(l — /i 2 )(ri ■ q r 2 • q) 2 . Averaging 
over the collisions that occur in a unit area per unit 
time gives the dissipation /. As the stress balances 
gravity, the work done by the second and third terms 
cancel, except for a contribution from a X yd x (vy} in the 
shear zones at the boundaries. This leaves the "heat" 
flow term (V • Q) to cancel the dissipation in the bulk. 
One must make a clear distinction here between "energy 
flux" and "heat flux" , something that has not been nec- 
essary in the diffuse regime JK|. From Eq. @), the 
energy flux from one grain to another in a collision is 



SF = i(l + /x)[(fi -q) 2 — (f 2 -q) 2 ]q. However, part of this 
energy flux is related to the coherent transfer of momen- 
tum (i.e. it is the work done by the stress tensor). As 
a result, the heat flux is the uncorrelated part of the en- 
ergy flux, Q = F + (a ■ v). With this done, one can relate 
the heat flux Q to the gradient of the granular temper- 
ature T = i((w 2 ) + ((v y - (v y )) 2 )) by Q = -kVT. As 
demonstrated in Figure ||(d), the conductivity k = pf c , 
not a terribly surprising result although different from 
what one hnds in the diffuse regime. 

In conclusion, we find that the dynamics of gravity 
driven flow within the regime studied do not fit into pre- 
vious descriptions of granular flow in the diffuse regime. 
While the model of instantaneous collisions with only 
normal forces is not necessarily general to all granular 
materials, much of the behavior, such as the presence of 
velocity fluctuations in the absence of shear does agree 
with experiments [Tq 1 ] on hard grains. While a complete 
theory requires more work, our results put fairly tight 
constraints on any such theory as well as providing very 
helpful information as to what are reasonable (and un- 
reasonable) assumptions in its formulation. 
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